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ABSTRACT 

Context. By means of chemical evolution models for galaxies of different morphological type, we have performed a detailed study of 
the evolution of the cosmic dust properties in different environments: the solar neighbourhood, elliptical galaxies and dwarf irregular 
galaxies. 

Aims. Starting from the same formalism as developed by Dwek (1998), and thanks to the uptodate observations available in the solar 
vicinity, we intend to study the effects of dust in the chemical evolution of different types of galaxies and at the same time to refine 
the investigation of the parameter space and achieving a satisfactory fine tuning of the parameters involved in our study. 
Methods. We have taken into account dust production from low and intermediate mass stars, supernovae II and la as well as dust 
destruction and dust accretion processes in a detailed model of chemical evolution for the solar vicinity. Then, by means of the same 
dust prescriptions but adopting different galactic models (different star formation histories and presence of galactic winds), we have 
extended our study to ellipticals and dwarf irregular galaxies. In all these systems, dust evolution has been calculated by means of 
chemical evolution models which relax the instantaneous recycling approximation and already reproduce the main features of the 
various galaxies. 

Results. We have investigated how the assumption of different star formation histories affects the dust production rates, the dust 
depletion, the dust accretion and destruction rates. We have predicted dust to gas and dust to metals ratios in very good agreement 
with those observed in the solar vicinity. We have shown how the inclusion of the dust treatment is helpful in solving the so-called 
Fe discrepancy, observed in the hot gaseous halos of local ellipticals, and in reproducing the chemical abundances observed in the 
Lyman Break Galaxies. Finally, our new models can be very useful in future detailed spectro-photometric studies of galaxies. 

Key words. ISM: dust, extinction; ISM: abundances; galaxies: abundances; galaxies: evolution. 



1. Introduction 

■ The presence of dust in local and distant galaxies is indicated by 
[ various observational evidences. The light emitted by stars inter- 
acts with the dust grains in a wavelength-dependent manner: the 
light emitted in the ultra-violet and optical bands is absorbed and 
scattered by the dust grains. This is the well-known phenomenon 
called dust extinction, which is taking place both inside and out- 
side galaxies (Aguirre 1999). In the optical band, the larger is the 
emission wavelength, the lower is the extinction effect caused by 
dust. 

In our Galaxy and in local starburst galaxies, dust grains are the 
main contributors to the emission in the mid and far infrared 
bands. In fact, the light absorbed by dust at UV and optical wave- 
lengths is then thermally re-emitted at much longer wavelengths, 
in the range 10-1000 fim. Also the chemical composition of the 
galactic interstellar medium (ISM) is strongly influenced by the 
presence of dust. Some chemical elements, called refractory, in 
the gas phase are subject to dust depletion and a fraction of their 
total abundance is incorporated into solid grains. Examples of 
refractory elements are Fe, Si, Mg, Ni. For these elements, the 
abundance measured in stellar surfaces is considerably higher 
than the gas phase abundances (Li 2005). 
Dwek (1998, hereinafter D98) has developed a chemical evolu- 
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tion model to study the dust content of the Milky Way galaxy 
and its evolution. He has focused on some refractory chemical 
species and for each one he has suggested a set of dust con- 
densation efficiencies, calculating the dust production rates from 
low and intermediate mass stars, type la and type II supernovae 
(SNe). For these elements, he has investigated also the dust de- 
struction and the accretion rates, the dust to gas ratios and the 
dust fractions at the present time. 

Since then, the amount of observational data regarding the 
local content of dust has noticeably grown (Draine 2003), allow- 
ing a reliable fine tuning of the parameters involved in theoretical 
dust evolution studies. By means of a one-zone chemical evo- 
lution for the solar neighbourhood, Zhukovska, Gail & Trieloff 
(2007) studied production by stars of various dust species. By 
considering the available observations of presolar dust grains 
in meteorites, their study has allowed them to put some con- 
straints on the dust condensation efficiencies from SNe and AGB 
stars. Moving to extragalactic objects, new observations call for 
a theoretical investigation of the dust content of galaxies of dif- 
ferent morphological types. In particular, the SCUBA camera 
has revealed a large number of dusty sources at high redshfit 
which might be young massive ellipticals undergoing a starburst 
(Lilly et al. 1999, Eales et al. 2000). For what concerns, in- 
stead, local ellipticals, new measurements have rejuvenated the 
interest in assessing the issue of the so-called iron-discrepancy 
(Arimoto et al. 1997) by means of the dust. In this paper, we 
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present new set of chemical evolution models, with updated nu- 
cleosynthesis prescriptions. All of these models are successful 
in reproducing the abundances of various galaxies (spirals, el- 
lipticals, irregulars). By means of our chemical evolution mod- 
els, we study dust evolution in various environments. Starting 
from the same formalism developed by D98, first we focus on 
the Solar Neighbourhood (S.N.) and we extend the analysis by 
D98, improving and deepening the investigation of the parame- 
ter space. We focus on the differential roles of type la and II SNe 
in dust destruction and production and the evolution of the ac- 
cretion and destruction rates. When possible, we compare these 
quantities with the data available in the literature computed by 
other authors. In this paper, we perform a chemical evolution 
study of various dust species in environments dilTerent than the 
Milky Way Galaxy. This is possible by taking into account in de- 
tail the stellar lifetimes, allowing us to predict the properties of 
dust species such as C and Fe, produced by stars with lifetimes 
spanning from ~ 0.03 Gyr up to ~ 10 Gyr. We focus on the ef- 
fects that different star formation histories have on the evolution 
of the dust content of galaxies of different morphological types. 
This paper is organized as follows: in section 2 we present the 
chemical evolution equations used to compute dust evolution 
and all the parameters involved. In section 3 we discuss our re- 
sults for the S.N., for elliptical galaxies and for dwarf irregular 
galaxies. Finally, in Section 4 we draw some conclusions. 



2. The chemical evolution model 

2.1. The chemical evolution equations 

By means of detailed chemical evolution models, it is possible 
to follow the evolution of the abundances of several chemical 
species and of the dust content of spirals, elUptical and irregular 
galaxies. In all models the instantaneous recycling approxima- 
tion is relaxed and the stellar lifetimes are taken into account. 
Detailed descriptions of the chemical evolution models can be 
found in Matteucci & Tornambe (1987) and Matteucci (1994) 
for elliptical galaxies, Chiappini et al. (1997, 2001) for the 
spirals and Bradamante et al. (1998) for irregular galaxies. 
In our picture, elliptical galaxies form as a result of the rapid 
collapse of a homogeneous sphere of primordial gas where 
intense star formation (SF) is taking place at the same time as 
the collapse proceeds. SF is assumed to halt as the energy of the 
ISM, heated by stellar winds and supernova (SN) explosions, 
exceeds the binding energy of the gas. At this time a galactic 
wind occurs, sweeping away almost all of the residual gas. After 
the SF has stopped, the galactic wind is maintained by type 
la SNe, and its duration depends on the balance between this 
heating source and the gas cooling (we refer the reader to Pipino 
et al. 2002, 2005). 

For spiral galaxies, the adopted model is calibrated in order 
to reproduce a large set of observational constraints for the 
Milky Way galaxy (Chiappini et al. 2001). The Galactic disc is 
approximated by several independent rings, 2 kpc wide, without 
exchange of matter between them. In our picture, spiral galaxies 
are assumed to form as a result of two main infall episodes. 
During the first episode, the halo and the thick disc are formed. 
During the second episode, a slower infall of external gas forms 
the thin disc with the gas accumulating faster in the inner than 
in the outer region ("inside-out" scenario, Matteucci & Francois 
1989). The process of disc formation is much longer than the 
halo and bulge formation, with time scales varying from ~ 2 
Gyr in the inner disc to ~ 7 Gyr in the solar region and up 
to 20 Gyr in the outer disc (see table 1). In this paper, we are 



interested in the study of dust evolution in the S.N. For this 
purpose, we focus on a ring located at 8 kpc from the Galactic 
centre, 2 kpc wide. 

Finally, irregular galaxies are assumed to assemble from infall of 
protogalactic small clouds of primordial chemical composition, 
until masses in the range ~ IO^Mq are accumulated, and to 
produce stars at a lower rate than spirals. 

All the models used in this paper consider only one gas 
phase. For a single-phase gas, the chemical evolution equation 
for a given chemical element i takes the following form: 



dGiit) 
dt 



= -mxiit) 

Jm, 



+A I <p(m) ■ 

-l-(l - A) ifrit - Tm)Qim{t - Tm)<l>{m)dm 
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where G,(0 = Mg(t)Xi(t)/M,o, is the gas mass in the form 
of an element ; normalized to a total fixed mass Mfot and 
G{t) = Mg(t)/M,o, is the total fractional mass of gas present in 
the galaxy at the time t. The same quantities can be defined in 
terms of the surface gas and mass densities, especially in spiral 
galaxies. X,(f) is defined as the abundance by mass (or mass 
fraction) of the element / (for a comprehensive discussion of 
this equation, see Matteucci & Greggio 1986). The quantity (//(t) 
is the star formation rate (SFR). The term -if>(t)Xi(t) gives the 
rate at which the element ; is subtracted from the ISM by the 
SF process. The second term is the rate at which each element 
is restored into the ISM by single stars with masses in the range 
Ml - Mb^„ , where is the minimum mass contributing, at a 
given time f, to chemical enrichment (the minimum is 0.8Mq) 
and Mb,„ is the minimum binary mass allowed for binary 
systems giving rise to type la SN (3Mq, Matteucci & Greggio 
1986). The quantities Qmi{t - t„) (where is the lifetime 
of a star of mass m) contain all the information about stellar 
nucleosynthesis for elements either produced or destroyed 
inside stars or both (Talbot and Arnett 1971). The third term 
represents the enrichment due to binaries which become type la 
SN, i.e. all the binary systems with total mass between Mb^ and 
Mbm - 16Mo. For the type la SN progenitor model, the Single 
Degenerate (SD) scenario is assumed, where a C-O white dwarf 
explodes by C-deflagration mechanism after having reached 
the Chandrasekhar mass (1.44Mo), owing to progressive mass 
accretion from a non-degenerate companion (Whelan & Iben 
1973). The parameter A represents the unknown fraction of 
binary stars giving rise to type la SN and is fixed by reproducing 
the observed present time SN la rate. In this third term, both 
quantities ^ and Qmi refer to the time t-tm^, where indicates 
the lifetime of the secondary star of the binary system, which 
regulates the explosion timescale. yu = M2IMB is the ratio 
between the mass of the secondary component M2 and the total 
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mass of the binary system Mb, whereas /(//) is the distribution 
function of this ratio. Statistical studies indicate that mass ratios 
close to 0.5 are preferred, so the formula: 



/(//) = 2i-^(i + r)/i^ 



(2) 



is commonly adopted, with j - 2 (Greggio & Renzini 1983). 
jUmm is the minimum mass fraction contributing to the SNIa rate 
at the time t, and is given by: 



f^min — fflQX 



(M2(t) M2- 0.5Mb 



I Mb 



Mb 



(3) 



The fourth term represents the enrichment due to stars in the 
mass range Mb^ - Mb„ which are either single, or, if in bi- 
naries, do not produce a SN la event. In this mass range, all 
the stars with masses m > 8Mq will explode as type II SNe, 
which in our picture are assumed to originate from core col- 
lapse of single massive stars. The fifth term represents the en- 
richment of stars more massive than Mb„, all of which explode 
as core collapse SNe. As the upper mass limits contributing 
to chemical enrichment, we assume Mu - lOOM©. Finally, 
the last two terms account for infall of external gas and for 
galactic winds, respectively. For the infaU term, an exponential 
law with different timescales is adopted for spirals and irreg- 
ulars (see Calura & Matteucci 2006b). Concerning ellipticals, 
Pipino & Matteucci (2004) showed that, in order to satisfy the 
largest number of photo-chemical properties, a quick infall, with 
a shorter timescale for the more massive objects, is needed. 
A non-zero outflow term is present in the equations describ- 
ing ellipticals and irregular galaxies. In both of these galaxies, 
a galactic wind develops as the thermal energy of the gas heated 
by SN explosions exceeds the binding energy of the gas (see 
Bradamante et al. 1998, Pipino & Matteucci 2004). The bind- 
ing energy of the gas is strongly influenced by assumptions con- 
cerning the presence and distribution of dark matter (Matteucci 
1992); for the model adopted here a diffuse iRe/Rd=0.l, where 
Re is the effective radius of the galaxy and Rd is the radius of 
the dark matter core) but massive {Mjark/Mium = 10) dark halo 
has been assumed (see Bertin et al. 1992). In the case of spiral 
galaxies, the outflow term is set to zero. 

2.1 .1 . The initial mass function 

(p(m) is the initial mass function (IMF), assumed to be constant 
in space and time and normalized to unity in the mass interval 
0.1 - IOOMq. For the spiral galaxies, we adopt a simplified two- 
slope approximation to the actual Scalo (1986) IMF, similarly to 
what is done in Matteucci & Fran9ois (1989), which is expressed 
by the formula: 



0Scalo('w) ■■ 



19 ■ m-'-35 
24 ■ m-'-™ 



if m < 2 Mq 
if m> 2Mq, 



whereas for elhpticals and irregulars we adopt a Salpeter 
(1995) IMF, of the form 

0Saip(m) = O.17-m-'-35 

The reason for such a choice reUes mainly on the chem- 
ical abundances and on the metal content observed in each 
morphological type (see Calura & Matteucci 2006b). In fact, 
it is well known that a Salpeter (1955) IMF can account for 



the abundances observed in local ellipticals and dwarf galaxies 
(Pipino & Matteucci 2004, Recchi et al. 2002), whereas it leads 
to an overestimation of the metal abundances in spiral discs 
(see Romano et al. 2005) which, on the other hand, are well 
accounted for by means of the Scalo (1986) IMF. 



2.1.2. Tine nucleosyntliesis prescriptions 

The yields used in the present work are separated into three 
groups: yields of low and intermediate mass stars, yields of type 

la SNe and yields of massive stars. 

Low and intermediate mass stars (i.e. with masses 
O.8M0 < m < 8M0) contribute through quiescent mass 
loss and planetary nebula phase to the ISM metal enrich- 
ment, type la SNe are assumed to originate from exploding 
white dwarfs in binary systems, characterized by total masses 
3Mq < nihin < 16Mq, according to the Matteucci & Recchi 
(2001) best model. Finally, we assume that single massive 
stars with initial masses in the range > 8 - lOOM© explode 
as core collapse SNe. The nucleosynthesis prescriptions are 
common to all models. For massive stars and type la SNe, we 
adopt the empirical yields suggested by Frangois et al. (2004), 
which are substantially based on the Woosley & Weaver (1995) 
and Iwamoto et al. (1999) yields, respectively, and are tuned 
to reproduce at best the abundances in the S.N. For low and 
intermediate mass stars, we adopt the prescriptions by van den 
Hoek & Groenewegen (1997). 



2.1 .3. Tine star formation rate 

For elhpticals and irregulars galaxies, the SFR ^(?) in our models 
is a Schmidt (1959) law expressed as: 



(4) 



with ^ = 1. Here Git) = X G,(0 is the normahzed total gas 
density or gas fraction at the time t. The quantity v is the effi- 
ciency of SF, namely the inverse of the typical time-scale for SF, 
and is expressed in Gyr~^. 

In the case of spiral galaxies, the SFR expression (Chiappini 
et al. 1997) is: 



(A(r,/) = v[ 



cr(r,t) 
o-(ro, t) 



-|2(yfc-l)|- 0"(^fa/) j<:-l^ 



o-(r,0 



ISM 



(r,t) 



(5) 



where v is the SF efficiency, cr(r, t) is the total mass (gas + stars) 
surface density at a radius r and time t, cr{rQ, f) is the total mass 
surface density in the solar region and crisMif, t) is the ISM sur- 
face mass density. For the gas density exponent k a value of 
1.5 has been assumed by Chiappini et al. (1997) in order to en- 
sure a good fit to the observational constraints for a large set 
of local spirals (Kennicutt 1998). The efficiency of SF is set to 
V = IGyr"', and becomes zero when the gas surface density 
drops below a certain critical threshold. For the SF, we adopt a 
threshold gas density CTth ~ IM^pc''^ in ffie disc as suggested by 
Kennicutt (1989). The difference between equations 4 and 5 is 
due to the different normalizations of the quantities involved in 
the chemical evolution equations (see section 2.1). Furthermore, 
the model used for spirals is a multi-zone one, where ffie SFR ex- 
pression is a function of ffie galactocentric radius. On the offier 
hand, the models of ellipticals and irregulars are one-zone and 
the SFR has a simpler expression. It is important to note that the 
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parameterization of the SFR in both equations 4 and 5 involves 
two different quantities, i.e. v and k, which act in the same way. 
In table 1 , we show the adopted parameters for all the chemical 
evolution models described in this section. In Figure[T] we show 
different SFRs for various models used in this work. In panel 
(a), we show the SFRs for the two models considered here for 
spiral galaxies: the S.N. model (solid line) and a model for the 
disk at a distance of 16 Kpc from the Galactic center (dashed 
line). Concerning the S. N. model, the most striking features are 
the SF hiatus at 1 Gyr (see Chiappini et al. 1997, 2001) and the 
threshold-dominated SF after 10 Gyr. Both features have impor- 
tant consequences on the behaviour of several quantities stud- 
ied in this paper. The model used for ellipticals is the one by 
Pipino et al. (2005). This is a first step in the self-consistent 
study of both optical and X-ray properties of elliptical galax- 
ies by means of a chemical evolution code. Detailed cooling 
and heating processes in the ISM are taken into account using 
a mono-phase one-zone treatment, allowing us a reliable mod- 
elling of the galactic wind regime. For reproducing a giant ellip- 
tical we make use of the Pipino et al. (2005) case Hal, which 
is characterized by a luminous mass Mj^m = lO'^M©, with a SF 
efficiency v = 25 Gyr ' and an infall timescale t - 0.2 Gyr, in 
which a mild secular gas accretion from the suiTounding IGM 
is allowed. For a more typical spheroid (L ~ L.) we refer to 
their model Lai, which has a luminous mass Mium - 10"Mo, 
a SF efficiency v = 15 Gyr ' and an infall timescale t - 0.3 
Gyr. We stress that, due to the uptodate nucleosynthesis pre- 
scriptions and to the inclusion of dust, the models presented here 
are intended as an improved version of the original Pipino et al. 
(2005) models Hal and Lai. To study the properties of a Lyman 
Break Galaxy, we use a model charactherized by a luminous 
mass Mium = lO'^M©, a star formation efficiency v = 5 Gyr"' 
and an infall timescale t = 0.5 Gyr, referred to as LBGal model. 
In our scheme, large galaxies form the bulk of their stars in a 
shorter timescale than small galaxies, according to the "inverse 
wind" scenario. As shown by Matteucci (1994), this scenario al- 
lows us to reproduce the coiTelation between the [Mg/Fe] ratio 
and the galactic mass observed in local ellipticals. 

In panel (b), we show the predicted SF history for two 
models of elliptical galaxies (Lai and Hal), characterized by 
very high SFR values (from ~ 200 to ~ 3000 M^lyr) and by a 
starburst lasting ~ 0.4 - 0.9 Gyr. After this strong event, the SF 
stops and a galactic wind sets in for several Gyrs, until when 
a diffuse and hot gaseous halo, which surrounds the galaxy, 
is formed (for a more detailed description of this model, see 
Pipino et al. 2005). 

Finally, in panel (c) of Fig. [T] we show two different models 
tested for an irregular galaxy. The solid line represents a dwarf 
irregular with continuous SF (IC), whereas the dotted line is a 
starburst irregular. For the IC model, we assume a continuous 
star formation with SF efficiency of v/c = 0.05 Gyr On the 
other hand, the starburst model is similar to the model B used by 
Lanfranchi & Matteucci (2003). The SFR of this model consists 
of 3 bursts with an efficiency vi,i,rsi - O-SGyr^^ . The times of 
occurrence of the bursts are 1 Gyr, 10 Gyr and 13 Gyr. The 
durations for each burst are 0.02 Gyr, 0.02 Gyr and 0.2 Gyr, 
respectively (for further details, see Lanfranchi & Matteucci 
2003). The IC and burst models have been used to reproduce the 
abundance patterns observed in Damped Lyman Alpha (DLA) 
systems (Calura et al. 2003, Dessauges-Zavadsky et al. 2004) 
and allow us to reproduce the present-day features of local 
dwarf irregulars and blue compact galaxies (Recchi et al. 2002, 
Lanfranchi & Matteucci 2003), respectively. This latter model 
should be regarded as representative of a blue compact dwarf 



galaxy. Our assumption is based on the results by Lanfranchi & 
Matteucci (2003), who have shown that the N/O, C/O, Si/O and 
O/Fe ratios observed in BCGs can be explained by a model with 
two to seven short bursts of SF with efficiencies in the range 
v = 0.1 -Q.9Gyr-^ 



2.2. Chemical evolution of tine Dust 

The chemical evolution of an element / in the dust is computed 
by using the formalism developed by D98. Let X^ustjiO be the 
abundance by mass of the element / in the dust and G(f) the ISM 
fraction at the time t, the quantity 



/(f) = Xdusui ■ Gif) 



(6) 



represents the normalized mass density of the element / at the 
time t in the dust. The time evolution of Gdusrjit) is calculated by 
means of the following equation: 



dGdustjit) 
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xSW/ 
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estr * accr 



dt 



(7) 



The main differences between eq.s [T]and |2]concern all the 
integrals in the right sides of both equations. These integrals, in 
the case of eq.|2] contain the quantities 5^.^ , 6'" and 6'' , i.e. the 
condensation efficiencies of the element / in stellar winds, type 
la and type II SNe (see section 2.2.1). These quantities represent 
the fractions of the element / which is condensed into dust and 
restored into the ISM by low and intermediate mass stars, type 
la and type II SNe, respectively. The seventh and eighth terms 
of eq. |2] represent the dust destruction and accretion rates, re- 
spectively. These terms depend on the quantities Tdestr and Tgccr, 
which represent the typical timescales for destruction and accre- 
tion, respectively. These two quantities are discussed in detail 
in sections 2.2.2 and 2.2.3, respectively. Finally, the last term of 
eq.|7]accounts for possible ejection of dust into the inter galactic 
medium (IGM) by means of galactic winds. This term is absent 
in the equation for the S.N. model, but is taken into account in 
the elliptical and irregular models. 

2.2.1. Dust condensation efficiencies 

We assume the dust condensation efficiencies suggested by D98. 
The choice of the following dust condensation efficiencies is 
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Fig. 1. Predicted SFRs as a function of time for different chemical evolution models. Panel (a): S.N. (solid line) and a model for the 

outer regions of the Milky Way disc (dashed line); panel (b): two different elliptical galaxy models, i.e. model Lai (solid line) and 
Hal (dash-dotted hue) of Pipino et al. (2005) ; panel (c): two different irregulars, one with continuous SF (sohd hne) and the other 
with 3 starbursts (dotted line). 



Parameters 


Model 


y 

(Gyr-') 


Tinf 

(Gyr) 


o-,„, 


k 


IMF 


Milky Way, S.N. 


1 


7 


54 


1.5 


Scalo 


Mill^ Way, 16 Kpc 


1 


20 


5 


1.5 


Scalo 


Model 


V 

(Gyr-') 


Tinf 

(Gyr) 


Mlum 
(Mo) 


k 


IMF 


Elliptical, Lai 


15 


0.3 


10" 


1 


Salpeter 


Elliptical, Hal 


25 


0.2 


10'^ 


1 


Salpeter 


Irregular, IC 


0.05 


10 


10" 


1 


Salpeter 


Irregular, burst 


0.5 


10 


10" 


1 


Salpeter 



Table 1. 

Adopted parameters for the galactic models used in this work. The various models are listed in column I. In columns 2, 3, 4, 5 and 6 we present 
the adopted parameters, i.e. the SF efficiency y, the infall timescale Ti„/, the total surface mass density cr,ot (in the case of the Milky Way) or the 
baryonic mass M(„„ (in the case of the elliptical/irregular model), the Schmidt law exponent k and the IMF, respectively. 



motivated by various lines of evidence, as discussed in detail 
in D98. In this paper, we focus on the main refractory ele- 
ments which are depleted into dust in the cold phase. These 
elements are C, O, Mg, Si, S, Ca, Fe. As an assumption, we 
consider that only these elements can be incorporated into dust 
grains. We assume that dust grains can be of two types: carbon 
(C) dust and silicate (Si) dust. We are aware that dust compo- 
sition can be more complicated. For instance, Mathis (1996) 
presented a model including also the possibihty of compos- 
ite grains containing both carbon, siUcates and oxides. Li & 
Greenberg (1997) presented a trimodal dust model with large 
sihcate core-organic refractory mantle dust particles, small car- 
bonaceous particles and Polycyclic aromatic hydrocarbon parti- 
cles, able to reproduce the interstellar extinction and polarization 
observational constraints. Furthermore, different dust structure 
would also have an impact on dust destruction and accretion (see 
Greenberg & Li 1999, Jones et al. 1996). However, the inclusion 
of more complicated dust types is beyond the aims of this paper, 
since it would increase the number of free parameters involved 
in our study. According to the formalism developed in D98, the 



contributors to the dust production are (a) low and intermedi- 
ate mass stars, (b) type la SNe and (c) type 11 SNe. We neglect 
the contributions from Wolf-Rayet stars and novae (D98), both 
of which are believed to represent unimportant sources of dust 
production. 

a. Low and intermediate mass stars - In these stars, dust is 
produced during the Asymptotic Giant Branch (AGB) phase 
(Ferrarotti & Gail 2006 and references therein). We assume 
that dust formation depends mainly on the composition of the 
stellar envelopes. If Xq and Xc represent the O and C mass 
fractions in the stellar envelopes, respectively, we assume that 
stars with Xq/Xc > 1 are producers of silicate dust, i.e. dust 
particles composed by O, Mg, Si, S, Ca, Fe. On the other 
hand, C rich stars, characterized hy Xg/Xc < 1, are producers 
of carbonaceous solids, i.e. carbon dust (Draine 1990). Being 
Mi^ej(m) and Mi dusAm) the total ejected mass and the dust mass 
formed by the stars as functions of the initial mass m for the 
element i, respectively, we assume that for stars with Xq/Xc < 1 

Mdus,,ci>n) = 4^ ■ [Mcejim) - 0.75Mo,,/m)] 
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with = 1 and 
MdustAm) = 0, 

for all the other elements. For stars with Xo/Xc > 1 in 
the envelope, we assume 

Mdust,c{m) = 



with 6^^ =1 for Mg, Si, S, Ca, Fe and 



with fii being the mass of the i element in atomic mass 
units. 

b. Type la Supernovae - for these systems, we assume 
Mdust.cim) = 6'^[Mej,cim)] 



with 6'^ = 0.5; 



MdustAm) = 6fMej,iim) 



with (5f = 0.8 for Mg, Si, S, Ca, Fe; 



Mdus,,o(m) = 16 Xi Mejjim)/ Hi 

b. Type II Supernovae - For type II SNe, we adopt the same 
prescriptions as for type la SNe. In a recent paper, Zhukovska et 
al. (2007) by means of a chemical evolution model for the solar 
neighbourhood study the available observations of presolar dust 
grains in meteorites. Their analysis points towards condensation 
efficiencies lower than the ones suggested by D98, in particular 
concerning the silicates. In this paper, we test various assump- 
tions concerning these parameters and study the effects also on 
systems other than the solar neighbourhood. 

2.2.2. Dust Destruction 

Dust destruction is primarily due to the propagation of SN shock 
waves in the warm/ionized interstellar medium (McKee 1989, 
Jones et al. 1994). Evidences for dust destruction come from the 
observations of high velocity clouds, where an anti-correlation 
between the depletion levels and the cloud velocities has been 
found (ShuU 1978, McKee 1989). This was interpreted as an ev- 
idence for grain destruction in SN shocks (D98). Following the 
suggestions by McKee (1989) and D98, for a given element ; the 
destruction timescale Tdestr can be expressed as: 



Tdestrj - (eMsm) 



R: 



' gas 
'SN 



(8) 



Hence the destruction timescale is independent from the dust 
mass. MsNR is the mass of the interstellar gas swept up by the 



SN remnant. For this quantity, McKee (1989) suggests a typi- 
cal value of MsNR ~ 68OOM0, which is in agreement with the 
results from numerical studies of SN evolution (Thornton et al. 
1998). As suggested by McKee (1989), typical values for the de- 
struction efficiency e in a three-phase medium as the present-day 
local ISM are around 0.2, hence we assume: 



eMsm = 0.2 x 68OOM0 = 1360Mo 



(9) 



RsN is the total SNe rate, including the contributions by both 
type la and type 11 SNe. No significant modifications are re- 
quired in a single-phase gas modelling (as we did in this paper, 
McKeel989). We note, however, that McKee (1989) estimates 
eMsNR ~ 70Mq for a hot and rarefied medium as the gaseous ha- 
los surrounding ellipticals, therefore we run two models in which 
eMsNR = I36OM0 until the galactic wind and then we have an 
instantaneous transition to eMsm ~ 70Mq. The pre- wind dust 
evolution is obviously unaffected, whereas it leads to substan- 
tial changes in the late stage of the galactic evolution. According 
to the galaxy mass, we call them Lal+MK and Hal+MK, re- 
spectively, and we will show that they might produce interesting 
results, even though they use an over-simplified treatment of the 
dust destruction in the hot gas. Finally, in order to render the 
cases Lai and Hal more realistic, we further modified the de- 
struction treatment by implementing a thermal sputtering term, 
which is thought to be the dominant source of dust destruction 
in hot plasmas. In particular, following Itoh (1989), we assume 
that, in a ~1 keV plasma, nearly 90% of the dust grains will 
evaporate by thermal sputtering in Tdestr^^j ~ 10^/ne(yr cm"-'), 
where the electron density n^. has been self-consistently evalu- 
ated at each timestep. This translates into a new destruction term, 
namely: 



Gdi 



'ust,i 



R 



Tde 



SN 



Gdi 



USt,l 



(10) 



According to the galaxy mass, we call the models featuring this 
particular term as Lal+Itoh and Hal+Itoh, respectively. For el- 
liptical galaxies, these two models will be regarded as the fidu- 
cial ones. According to Itoh (1989), this term dominates the Fe 
grains evolution at late times. 

In section 3.1.2, we test the parameter e and we tune it in 
order to reproduce the depletion pattern observed in the S.N. 
Further refinements are presented in section 3.2. 

2.2.3. Dust Accretion 

Dust accretion occurs in dense molecular clouds, where volatile 
elements can condensate onto pre-existing grain cores, originat- 
ing a volatile part called mantle (D98, Inoue 2003). Direct ev- 
idences for dust accretion come from the observed large varia- 
tions of the depletion levels as a function of the density (Savage 
& Sembach 1996) and from the observed infra red emission of 
cold molecular clouds (Flagey et al. 2006), which is character- 
ized by the absence of small grain emission. These features can 
be accounted for by the coagulation of small grains on and into 
larger particles. Indirect evidence for dust accretion comes from 
the estimation of the grain lifetimes, which would be very small 
if no process could allow the grains to recondense and grow 
(McKee 1989, Draine & Salpeter 1979). For a given element 
i, the accretion timescale Taccr can be expressed as: 



Taccr = To,,7(l - fi) 

where 

Gdust,i 



fi 



Gi 



(11) 
(12) 
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According to eq. (TT] the accretion timescale is an increasing 
function of the dust mass. For the timescale tqj, typical values 
span from ~ 5 x 10^ yr, of the order of the lifetime of a typical 
molecular cloud (D98), up to ~ 2 x 10^ yr (D98). In this paper, 
we assume that the timescale tq,/ is constant for all elements, 
with a value of 5 x 10^ yr. 

In elliptical galaxies, we assume that dust accretion occurs only 
during the starburst epoch, when large amounts of cold gas and 
molecular H are available. After the onset of the wind and the 
end of the starburst, we assume that no molecular gas is present, 
hence no more accretion can occur. 

The observed molecular H content in dwarf irregular galaxies is 
very small, with molecular-to-atomic gas fractions of ~ 10% or 
lower (Lisenfeld & Ferrara 1998, Clayton et al. 1996). Motivated 
by these observational results, we assume that no accretion can 
occur in iiTegular galaxies. 

3. Results 

3. 1 . Dust evolution in the Solar Neighbourhood 

3.1.1. Dust Production Rates 

In Figure|2l we show the predicted evolution of the carbon (lower 
panel) and silicate (upper panel) dust production rates calcu- 
lated by means of the chemical evolution model for the S.N. 
In particular, we show the contributions to the rates by various 
sources: the low and intermediate mass stars (LIMS, solid lines), 
the type II SNe (dotted lines) and the type la SNe (dashed lines). 
Concerning the C dust, its production is dominated by LIMS 
throughout most of the cosmic time. At the present day (To ~ 13 
Gyr) a significant contribution is coming also by type II SNe, 
with type la SNe playing a negligible role. The frequent discon- 
tinuities in the lines are due to the effect of the star formation 
threshold. 

The production of silicate dust is dominated by type II SNe 
throughout most of the time. At late times, the contributions by 
type la SNe and type II SNe are comparable. On the other hand, 
LIMS are negligible contributors to the Si dust production rate. 

3.1 .2. The dust fractions 

An important test for the robustness of our formalism is repre- 
sented by the study of the metal fractions in dust. These quanti- 
ties are the ratios between the amount of a given element locked 
into dust and its total abundance and can be expressed as: 

^ ^ (13) 

J. 

Of particular interest is the comparison between the pre- 
dicted dust fractions for various elements and the observed ones. 
This analysis is very useful to test the various parameters in- 
volved in our study, in particular the dust condensation efficien- 
cies S^^,6''', and 6'', the dust destruction efficiency e and the 
dust accretion timescale to,,-, as well as the sensitivity of the 
results on the assumptions about these parameters. In Fig. [3] 
we show the predicted present-day fractions in dust for the el- 
ements studied in this work, compared to the values observed 
by Kimura, Mann & Jesseberg (2003) in the Local Interstellar 
Cloud. It is worth to note that the Local Interstellar Cloud con- 
sists of a warm medium, characterized by a temperature of 
~ 60007*^, whereas our models provide the depletion of the cold 
gas, with temperatures lower than ~ IQQK. In this work, the dust 



fractions observed by Kimura et al. (2003) are taken as refer- 
ence. From the logarithmic depletions 6 plotted in Figure 2 of 
Kimura et al. (2003), it is possible to derive the dust fractions of 
the cold gas through the formula 

/ = 1 - 10* (14) 

For C, O and S, the logarithmic depletions of the cold medium 
are very similar to the ones of the Local Intersetllar Cloud. For 
Si, Mg and Fe, the dust fractions of the cold medium are higher 
than the ones of the warm medium and are within the error bars 
plotted in Figure [3] In the two panels of Fig. [3] the predictions 
are obtained by testing two different assumptions concerning the 
dust condensation efficiencies. In the lower panel of Fig. [3] we 
present the predicted values calculated considering only dust 
production in stars. The open stars are the fractions calculated 
by adopting for the dust condensation efficiencies the prescrip- 
tions suggested by D98. The open pentagons are the fractions 
calculated by assuming for the dust condensation efficiencies 
a constant value of 0.1, regardless of the element and of the 
source of production. Following the prescriptions suggested by 
D98, a higher amount of metals is locked into dust than with 
the other choice. The only source of dust destruction present in 
this case, i.e. astration, does not seem to play a dominant role 
in determining the dust depletion pattern, once destruction by 
SNe and accretion in the ISM are neglected. We outline that 
neglecting these two processes is unrealistic since, as stressed 
in sections 2.2.2 and 2.2.3, there are robust evidences that de- 
struction and accretion play a non negligible role in the evolu- 
tion of dust in the S.N. (see McKee 1989, Draine 1990, Tielens 
1998). For purpose of comparison, we show also the observa- 
tional values derived by Kimura et al. (2003) (open squares with 
solid and dashed error bars). The chemical composition of dust 
grains cannot be directly observed. The only way to estimate 
the abundance of an element locked up into dust is by means 
of its gas phase abundance, which is directly observable, and 
by means of an appropriate value representing the cosmic abun- 
dance. Observationally, for a given element /, the abundance in 
dust is derived by means of the subtraction between the cosmic 
abundance (Xj c)obs and the observed gas abundances (Xgas,dobs- 

(X-i^dust)obs — {Xi^c)obs ~ (Xgas.dobs (15) 

The observational determination of the dust fractions de- 
pends on the assumptions on the total cosmic abundances 
iXj^c)obs and on the H ionization fractions xh = 0.25. In the 
two panels of Fig. |3] we consider four different sets of observed 
dust fractions (for details see caption of Fig. |3]l, calculated with 
various assumptions of both (X, c)oij ^^nd^H - 0.25. 
In the upper panel of Fig. [3] the predicted values have been cal- 
culated by taking into account dust production in stars, dust de- 
struction and dust accretion in the ISM. The calculation of dust 
destruction and accretion depends on the choice of two parame- 
ters, i.e. the destruction efficiency e and the accretion timescale 
To,,, respectively. In this case, the predictions have been calcu- 
lated by assuming e ~ 0.2 and To,,- = 5 x 10^ yr for all elements. 
For this choice, the dust fractions are nearly independent from 
the choice of the dust condensation efficiencies. This confirms 
the results by Dwek (1998), who showed that, at the present 
time, the accretion rate balances the destruction rate, and that 
the balance determines the depletion. Similar conclusions are 
drawn from the analysis of Zhukovska et al. (2007). By means 
of a different approach, Tielens (1998) attempted to determine 
the accretion and destruction rates from the observed depletions. 
The implication of these results is that, to derive constraints on 
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Fig. 2. Predicted dust production rates from various sources for a chemical evolution model of the S.N. In the lower (upper) panel, 
we show the results for the carbon (silicate) dust. Solid lines: contribution by low and intermediate mass stars (LIMS). Dotted lines: 
contribution by type II SNe. Dashed lines: contribution by type la SNe. 



the condensation efficiencies, one has to study the depletion pat- 
tern in systems where either destruction or accretion are absent. 
In section 3.3.4, we will show that systems suited to this study 
are the dwarf irregular galaxies, where dust accretion is likely to 
play a negligible role. 

From the upper panel of Fig. [3] we note also that for some ele- 
ments, the assumption of a constant destruction efficiency pro- 
vides a very poor fit to the observed dust fractions. These ele- 
ments are C, O and S, for which the predictions indicate present- 
time dust fractions of ~ 0.9, higher than the values observed in 
the Local Interstellar Cloud. 

The observed dust fractions can be reproduced in two ways, 
by assuming that either the dust accretion timescale or the the 
destruction efficiency depends on the physical properties of the 
chemical element. The dust accretion timescale is related to 
the micro-physics of the molecular clouds, which can not be 
modelled properly with our instruments. On the other hand, the 
destruction efficiency is connected to the heating of the ISM, 
mainly due to SNe and their explosion rate, which is accounted 
for in detail by our models. In this paper, we assume that the dust 
destruction efficiency depends on the properties of the chemical 
element. A physical justification is related to the condensation 
temperatures of the various elements. These quantities express 
the volatility of the elements in dust (Lodders 2003 and refer- 
ences therein). In general, elements with higher condensation 
temperatures are more likely to aggregate into dust grains and 
are more resistant to destruction. Lodders (2003) has calculated 
the condensation temperatures for various elements for a solar- 
system composition gas. Her results indicate condensation tem- 
peratures of 78 K, 182 K and 704 K for C, O and S, respectively. 
On the other hand, the condensation temperatures calculated for 
Fe, Si and Mg are 1357 K, 1529 and 1397, respectively, hence 
considerably higher than for C, O and S. Here, we assume that 
the destruction efficiencies vary as a function of the condensa- 
tion temperatures of the elements. For C, O and S we assume 
a dust destruction efficiency of ec,o,s - 0.8, whereas for Fe, Si 
and Mg we assume eFe.Mg,si - 0.2. As can be seen from Figure 
m where we show the dust fractions calculated by assuming the 
above set of dust destruction efficiencies, with these values it is 



possible to reproduce the observed dust fractions with good ac- 
curacy. 

In Figure |5] we show the time evolution of the dust fractions 
for various elements (C, O, Fe and S). We have omitted the Mg 
and Si because for both elements, the evolution of the fraction 
in dust is identical to the one predicted for Fe. The curves for C, 
O and S show a similar behaviour after the minimum at 2 Gyr, 
with an increase up to the present time which is steeper than the 
one for Fe. This fact is due to the different destruction efficien- 
cies adopted for the C, O and S group and for the Si, Fe and Mg 
group. 

3.L3. The dust destruction and accretion rates 

Once we have tuned the destruction efficiency e and chosen a 
realistic value of the accretion timescale tqj, we study the evo- 
lution of the accretion and destruction rates and we compare our 
values to the ones estimated by other authors. In the left panel 
of figure |6l we show the predicted evolution of the dust destruc- 
tion and accretion rates for the S.N. The evolution of both rates 
is sensitive to the SF and gas accretion history of the model for 
the S.N. In the literature, no observational estimate can be found 
for the accretion or destruction rates in the S.N. These quan- 
tities have been calculated by various authors on the basis of 
theoretical investigations (McKee 1989, Draine 1990, Jones et 
al. 1994, Tielens 1998 ). Our model is one-phase and represents 
one single medium dominated by cold, neutral gas. The rates 
predicted by our models with the average estimates by Draine 
(1990) and Tielens (1998) for a cold neutral medium. The values 
by Draine (1990) and Tielens (1998) are order of magnitude esti- 
mations, and are provided with no error bar. For the dust destruc- 
tion rate, we predict present-day values of ~ 5 x 10^'^yr"' and 
~ 8 X lO^^yr"' for silicate and carbon dust, respectively. Similar 
values are predicted for the accretion rates for both dust types. 
Our values are in good agreement with the order-of-magnitude 
estimations determined by Draine (1990) and Tielens (1998) 
with different approaches. 

Several recent observational and theoretical results have ad- 
dressed the role of SNe in dust production in the Milky Way 
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Fig. 3. Fractions in dust for various elements. Open stars: predicted present-day fractions calculated by adopting for the dust con- 
densation efficiencies the prescriptions suggested by D98. Open pentagons: predicted present-day fractions calculated by assuming 
for the dust condensation efficiencies a constant value of 0. 1 . The solid and open squares are the fractions observed by Kimura et al 
(2003) in the Local Interstellar Cloud using the set of cosmic abundances specified in their table 2 and 3, respectively. The squares 
with solid and dashed eiTor bars have been calculated assuming a H ionization fraction of xh = 0.25 andxH = 0.45, respectively. 
In the lower panel, the predicted values are calculated considering only dust production in stars. In the upper panel, the predicted 
values are calculated considering dust production in stars, dust destruction by SNe and dust accretion in the ISM. 
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Fig. 4. Fractions in dust for various elements. The solid pentagons are the predicted present-day fractions calculated by assuming a 
dust destruction efficiency of e = 0.2 for all elements. The open diamonds are the present-day fractions calculated assuming a dust 
destruction efficiency of ec,o.5 = 0.8 for C, O and S and whereas for efe,Mg,si - 0.2 for Fe, Si and Mg. The squares with solid and 
dashed error bars are the same as in Figure |3] 



(Dunne et al. (2003) and in high redshift galaxies (Hughes et al. 
1998, Bertoldi et al. 2003, Morgan & Edmunds 2003, Maiolino 
et al. 2004). Beside producing dust, SNe are the main sources of 
energy injection into the ISM and their hot cavities are known to 
destroy dust grains (D98). Our aim here is to compare the type 
la and II SN destruction and production rates and see how do 
these quantities evolve with cosmic time. In the right panels of 
Figure |6l we study the balance between production and destruc- 
tion by SNe for C dust (lower panel) and Si dust (upper panel). 



The thin dotted and solid lines are the predicted evolution of the 
dust production rates by type la and type II SNe, respectively. 
The thick dotted and solid lines are the predicted evolution of the 
dust destruction rates by type la and type II SNe, respectively. 
For type II SNe, dust destruction dominates over dust produc- 
tion throughout almost all the cosmic history. Only for a very 
short time (~ 0.3 Gyr after the beginning of SF), the type II SN 
production rate was larger than the type II SN destruction rate. 
This is in agreement with the results by D98, who found that at 
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Fig. 5. Time evolution of the fractions in dust for various elements: C (solid line), O (dotted line), Fe (dashed line) and S (dot-dashed 
line). The evolution of the dust fractions for Mg and Si is identical to the Fe one. 




lime (Gyr) 

Fig. 6. Left: Evolution of the dust accretion (upper panel) and destruction (lower panel) rates in the S.N. The solid and dotted lines 
represent the calculated rates for carbon and silicate dust, respectively. The open triangles and solid triangles are the present-day 
average values estimated by Draine (1990) and Tielens (1998), respectively, for a cold neutral medium. Right: Evolution of the 
production and destruction rates by SNe for C dust (lower panel) and Si dust (upper panel). The thin dotted and solid lines are 
the predicted dust production rates by type la and type II SNe, respectively. The thick dotted and solid lines are the predicted dust 
destruction rates by type la and type II SNe, respectively. 



earlier epochs, when the ISM metallicity was below Z ~ 0.001 
(reached, according to our model, in ~ 0.35 Gyr) SNe inject 
more dust into the ISM than they can destroy. 
On the other hand, for type la SNe during the first Gyr of galactic 
evolution. Si dust production has dominated over dust destruc- 
tion, but with negligible consequences on the total Si dust mass. 
We conclude that, according to our results, throughout the whole 
history of the S.N., the presence of SNe, in particular of type II 
SNe, is fundamental for the process of cycling the dust back into 
the ISM. 



3.1 .4. Dust to gas and Dust to metals ratios and their 
evolution 

We define the dust to gas ratio D as the ratio between the dust 
mass and the total ISM mass surface densities. 



D : 



O" dust 
O" ISM 



(16) 



while the dust to metal ratio Dz is the ratio between the dust 
mass and the metal mass surface density in the ISM, 



O" diisl 



D 



(17) 
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Fig. 7. Upper panel: predicted evolution of the dust to gas ratio in the S.N. (solid line) and at 16 Kpc from the Galactic centre 
(dashed line). The solid square with error bar is the observational value for the Galaxy taken from Issa et al. (1990). Lower panel: 
predicted evolution of the dust to metals ratio in the S.N. (solid line) and at 16 Kpc from the Galactic centre (dashed line). The open 
square with error bar is the observational value for the Galaxy assuming the dust to gas ratio estimated by Issa et al. (1990) and a 
solar metallicity of Zq = 0.0133 (Lodders 2003). 



Issa, MacLaren and Wolfendale (1990), by studying the D 
quantity in local spirals, have found a correlation between D and 
the metallicity. The implication of this fact is that the dust to met- 
als ratio is nearly constant for local spirals. In the upper panel of 
Figure |2l we show the predicted time evolution of D for two 
different zones of our Galaxy: the S.N. (solid line) and a cir- 
cular ring 1 kpc wide located at a distance of 16 Kpc from the 
Galactic center (dashed line). The latter represents a model for 
the outermost regions of the spiral disc (see table 1). The com- 
parison between these two models is interesting to understand 
the differences of the dust properties for two regions located in 
different zones of the Galactic disc. We plot also the Galactic 
value (solid square with error bars) Dq = 0.008, taken from Issa 
et al. (1990), which is well reproduced by our model. 

In the lower panel of Figure |2l we plot the time evolution of 
the dust to metals ratio for the S.N. (solid line) and the spiral 
outskirts (dashed line) models. We compare the predictions to 
an observational value derived for the Galaxy. By assuming for 
the S.N. a metallicity of Zq = 0.0133 (Lodders 2003), for our 
Galaxy we obtain Dz.q = 0.6. Interestingly, this value coincides 
with our predicted dust fraction for O (see Fig. |5]l, which is the 
element dominating the total metallicity. 

Chemical evolution studies have shown that the dust to met- 
als ratio should not vary much during galaxy evolution (Inoue 
2003, Edmunds 2001). Our predictions indicate that in the initial 
phases of galaxy evolution, in particular from the halo phase to 
the SF hiatus, Dr experiences strong variations, evolving rapidly 
from 0.15 to 0.8. After the peak at 1 Gyr, drops sharply and 
its evolution in the disc phase is much slower, i.e. of a factor of 
~ 2 from 2 Gyr up to the present time. This is consistent with 
the results by Inoue (2003). The present-day value is nicely 
reproduced by our models. The spiral outskirt model presents 
D- values higher than the S. N. model. Dust accretion is likely 
to take place in the outermost regions of the Milky Way owing 
to the presence of molecular clouds. In fact, recent near Infrared 
observations have revealed the presence of molecular clouds in 
regions located at galactocentric distances out to R > 13.5-20 



kpc (Kobayashi & Tokunaga 2000, Snell et al. 2002, Nakagawa 
et al. 2005). 

In table 2, we present a detailed comparison between the nu- 
meric values as predicted by means of our model for the S.N. 
with the predictions by D98, as well as with the available cor- 
responding observational results. Concerning the predicted ISM 
properties, our results are very similar to the ones of D98 and 
well within the available observational values. The small differ- 
ences are due mainly to different prescriptions for the SF law 
(see our equation 5 and eq. 15 of D98) and for the infall law, for 
which we assume the two infall expression by Chiappini et al. 
(1997), whereas D98 assumes a simple exponential law. 
Differences in the C and Si dust production rates from LIMS and 
SNe (la-HlI) are mainly due to the different nucleosynthetic pre- 
scriptions used in this paper and by D98. In particular, for stars 
of low and intermediate mass we use the metallicity-dependent 
yields by Van den Hoeck & Groenewegen (1997), whereas, for 
stars in the same mass range, D98 used the yields computed 
by Renzini & Voli (1981). For the dust production rates, in ta- 
ble 2 we decided not to show any observational value since the 
only existing observational constraints are the ones by Jones & 
Tielens (1994), which are provided without any error bar. 
Differences in the dust accretion rates DAR are due to the differ- 
ent present day gas densities and to the different nucleosynthesis 
prescriptions, whereas the small differences in the destruction 
rates DDR are due to the different star formation laws, which 
cause differences in the predicted SNRs. Finally, in the last two 
lines of table 2, we compare the C and Si dust to gas ratios as 
computed by us at the time when the solar system formed (i.e. at 
9.5 Gyr), with the values as predicted by D98. 



3.2. Dust evolution in elliptical galaxies 

In this Section we make use of the best values for the parameters 
derived above for the S.N., in order to extend our formalism to 
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Table 2. Comparison between some present-day output quantities (listed in the first column) predicted by means of our model for 
the solar neighbourhood (second column), with the same quanties as computed by D98 (third column) and with the corresponding 
observational values (when available, fourth column). Rq - %kpc is the solar distance from the Galactic center. References; a: Rana 
1991. b: Dickey 1993. c: Kulkarni & Heiles 1987. d: Portinari et al. 1998. 



Quantity 


Our model 


D98 model 


Observed 


ISM properties 






2.7 


3.8 


2-10" 


7 


8.8 


7-16 




1 


0.8 


0.3-1.5'' 


Tiiict ni"r\rinr'ti r\ri in ctncc 
LvLlSL UHJU.U.L.LHJ1I 111 SLdlS 




DPRcLiMS (Rq, To){MQpc-^Gyr-^) 
DPRc,sn(Ro, To)(Mapc-^Gyr-') 


4 X 10-* 
6 X 10-* 


2.8 X 10-3 
1.6 X 10-3 




DPRs.MMS (Re, To)(Mepc-^Gyr-') 


1.4 X 10-3 


3.7 X 10-3 




DPRs,.sN(Re, To)(Mepc-^Gyr-') 


8 X 10-3 


1.1 X 10-2 




Dust accretion in the ISM 




DARciRe, To){Mepc-'Gyr-') 


0.065 


0.035 




DARs.iRe, To)(Mepc-^Gyr-') 


0.23 


0.18 




Dust destruction due to SN shocks 




DDRciRe, To)(Mepc-'Gyr-') 
DDRs,(Re,To)(Mepc-'-Gyr-') 


0.058 


0.033 




0.25 


0.19 




Dc 


8.5 X 10-* 


1.7 X 10-3 






4.8 X 10-3 


7 X 10-3 





the class of early-type galaxies, implementing eq.|2]in the model 
by Pipino et al. (2005). 

As a first step, we do not modify the parameters related to 
processes such as the dust growth, which are mainly governed by 
the microphysics, and thus they should not depend on the galac- 
tic morphology. However, after the onset of the galactic wind, el- 
lipticals are basically devoid of cold gas. Therefore, we assume 
that, during this phase, no dust accretion can occur in elliptical 
galaxies. 

Concerning the last term in eq. [T] now dust is allowed to 
escape the ISM during the galactic wind, although no conclusion 
can be drawn on its fate. Therefore, the dust loss rate in eq. [T] 
switches on during the galactic wind and it has been assumed to 
follow the gas flow rate. In particular, we recall that Pipino et 
al. (2002, 2005) successfully reproduced the amount of metals 
injected into the intra cluster medium (ICM) by assuming that, 
at each timestep, the ejected mass A M was set by the following 
condition: = where AE was the difference between the 

thermal energy E^h and the gas binding energy. Here we assume 
a similar condition to hold also for dust flow rate. Therefore, we 
have that the fraction of the ISM (i.e. dust plus gas) which can 
escape, scales with AE as: 

AM _ AMd„„ + AMg„s _ AE 

We are aware that the infrared emission from grains might 
change the cooling function (Draine 1981). To tackle this is- 
sue, it is necessary to study interstellar dust from a spectro- 
photometric point of view. This will be the subject of a forth- 
coming paper. 

3.2.1. The dust production rates 

In Figure [8] we show the predicted dust production rates for 
the Lai model, describing an elliptical galaxy of stellar mass 
~ 1O"M0. During the starburst phase, lasting ~ 1 Gyr in this 
case, the major C and Si dust producers are type II SNe. After 



the end of the starburst. Si dust production is dominated by type 
la SNe, whereas C dust production is mainly due to type la SNe. 
At late times, i.e. after 10 Gyr, LIMS are the major C dust pro- 
ducers. 

3.2.2. The evolution of the dust fractions 

In Figure |9] we show the predicted evolution of the dust frac- 
tions locked up in C, O, S and Fe for the fiducial case Lal+Itoh. 
In the very early phase of the evolution (i.e. t <200 Myr, when 
the gas is still infalling and the SFR has just reached its max- 
imum intensity), we predict that nearly 80% of Si and Fe are 
locked into dust, whereas only 40% of C and 10% of O. This 
finding is in rough agreement with the depletion pattern found 
for spirals. After an initial peak, the dust fractions undergo a 
sudden decrease due to the high SN activity, which leads to high 
destruction rates. The decrease in the dust fractions is then en- 
hanced by the mass removal as soon as the galactic wind devel- 
ops. After the end of the galactic wind (around 7 Gyr), the dust 
mass can increase again, before the sputtering due to the hot X- 
ray emitting gas can have a major role in the grain destruction. 
An interesting case is represented by the S dust fraction. In spi- 
rals, the evolution of the S dust fraction is similar to the ones 
of C and O, mainly because of the similar destruction efficien- 
cies. In these systems, the destruction is very efficient because 
of the continuous explosions of type II SNe. On the other hand, 
dust destruction plays a minor role in ellipticals, presenting long 
periods of no star formation activity. In these systems, the evo- 
lution of the S dust fraction is close to the one of Fe, since the 
two elements have similar condensation efficiencies. 

3.2.3. The dust fractions in present-day ellipticals: local 
ellipticals and the Fe discrepancy 

The main answer that we expect from our analysis is whether the 
inclusion of the dust treatment may help in solving the so-called 
iron discrepancy (Arimoto et al. 1997) in X-ray spectra of the hot 
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halos surrounding ellipticals. When the ASCA satellite provided 
the first reliable measure of the iron abundance in the hot ISM 
of ellipticals (e.g. Awaki et al. 1994; Matsumoto et al. 1997), 
it was much lower than the solar value, at odds not only with 
theoretical models for elliptical galaxies available at that time 
(Arimoto & Yoshii 1987; Matteucci & Tornambe 1987), which 
predicted that their ISM should exhibit [Fe/H]> 0, but also with 
the mean metallicity of the stellar component inferred from op- 
tical spectra. This issue then had been partly alleviated by tak- 
ing into account temperature gradients (Buote & Fabian 1998; 
Buote 1999). Recently, Humphrey & Buote (2006) reported 
emission-weighted Fe abundances up to ~ 2-3 solar in a sample 
of 28 early-type galaxies, corresponding to [Fe/H]=0.3 - 0.48. 
Nevertheless, the iron discrepancy persists when comparing the 
recent results by Pipino et al. (2005), predicting [Fe/H]> 0.85 
in the ISM of elliptical galaxies, to the most recent abundance 
measurements in the X-ray spectra. 

Fe condensation in dust is among the possible physical 
mechanisms often invoked to solve this issue (see Arimoto et 
al. 1997 for a comprehensive analysis). Recent far-infrared ob- 
servations, in fact, claim that the dust mass in ellipticals could 
be ~ IO^^^Mq, higher than previous estimates (Temi et al. 2004) 
by a factor of ten, with ~ 2 ■ lO^M© of dust residing in the very 
central galactic regions (Leeuw et al. 2004). Independent recent 
observations carried out by means of the SCUBA camera on lo- 
cal ellipticals (Vlahakis et al. 2005) seem to confirm these val- 
ues, although the sample should be enlarged and the effects of 
synchrotron radiation clarified. Recently, several observationalal 
determinations of the dust masses in ellipticals have been possi- 
ble thanks to the Spitzer Space Telescope. 
In elliptical galaxies, Kaneda et al. (2007) derive dust masses 
typically of the order of 10^ - 10^ Mq. In an early type galaxy 
of the Virgo group, Panuzzo et al. (2007) obtain a dust mass 
of 1.2 X IO^Mq. These new observational values are well re- 
produced by our fiducial models, which predict a dust mass 
of ~ 1.1 ■ 10'' for more than 10 Gyr old giant ellipticals 
(~ 0.8 ■ IQ^Mq of dust for a more typical galaxy like our model 
Lai). 

The predicted evolution of the Fe dust destruction rates and of 
the Fe fraction in dust for the Lai model assuming the Itoh and 
the McKee prescriptions for the destruction are presented in the 
upper and lower panel of Figure [TOl respectively. In Figure [TT] 
we show the predicted evolution of the Fe abundance in the hot 
phase of a typical elliptical galaxy by adopting the Itoh (1989, 
solid line) and the Mckee (1989, dotted line) prescriptions. For 
the fiducial model, the predicted present-day iron abundance in 
the hot gas is [Fe/H]=-i-0.6 dex (for both Lai and Hal). These 
abundances are a factor of 2-5 less than the values predicted 
by the original Pipino et al. (2005) models Lai and Hal, re- 
spectively. Several factors act in determining this result: dust 
depletion (up to 20% of Fe is in dust at late times), a more 
prolonged wind phase with respect to the original model Hal 
(which leads to a more efficient Fe ejection), different stellar 
yields. Moreover, since the Fe emission line dominates the X- 
ray spectrum of the hot ISM, a reduced Fe abundance in the gas 
leads to a better agreement between the predicted X-ray lumi- 
nosity (2.1 ■ 10'*^erg sec"') and temperature (~ 0.7 kev) for the 
hot halo of giant ellipticals with the observed values, than in 
Pipino et al. (2005) case. When adopting the Lai -hM/T prescrip- 
tion (dotted lines in Figs. 1121 and [TOl l. the grain destruction is 
considerably reduced and up to ~ 4 ■ IQ^Mq of dust can survive 
after 12 Gyr of galactic evolution. At this stage, the amount of 
Fe locked up in dust is close to 60% (for a comparison, the frac- 
tion for O and C are only 25% and 12%, respectively), leading 



to a [Fe/H]=0.15 (Hal) - 0.3 (Lai) dex as well as [O/Fe]=-0.12 
dex, in excellent agreement with the available data (see Pipino 
et al. (2005), table 4 for a compilation of the most recent ones): 
this might be the solution of the iron discrepancy. Contrary to 
other authors (Arimoto et al. 1997), we suggest that the iron dis- 
crepancy might be solved by taking into account the dust con- 
densation of Fe. New and elliptical-dedicated calculations of the 
sputtering by hot (i.e. T > lO^K) plasmas are needed to eventu- 
ally asses this issue. 

3.2.4. The dust flow into the IGM/ICM 

The amount of dust ejected by a giant elliptical galaxy can be ob- 
tained by integrating eq[T8]over time. The amount is ~ 3 ■ 10^ Mq 
for the fiducial Hal model. For the model Lai (fiducial case), in- 
stead, the ejected dust mass is ~ 5 ■ 10*^ Mq. Once the dust grains 
are ejected into the IGM, no firm conclusions on their conse- 
quent evolution can be drawn. For instance, we need a detailed 
comparison of the efficiency of dust removal by means of the 
SN-driven wind with the dust destruction rate due to the SN ex- 
plosions. Moreover, a careful treatment of the dust survivability 
in a medium such as the hot ICM/IGM is beyond the scope of 
the paper 

3.2.5. The dust fractions in young ellipticals: hints from 
Lyman Break and SCUBA galaxies 

Unfortunately, we do not have direct measurement of the dust 
composition in young ellipticals. Therefore, in this section, we 
will try to constrain our models by means of the objects which 
are believed to be their high-redshift counterparts. 

It has been argued that the so-called SCUBA sources may be 
ellipticals in the process of formation (Lilly et al. 1999; Eales et 
al. 2000). The inferred SFRs of ~ 100 - lOOOMo/yr are compa- 
rable to those assumed by our models, as well as the estimated 
duration of the formation process and stellar mass involved. A 
useful measure of the dust mass in such high redshift objects is 
given by Dunne et al. (2003). They fitted the dust mass function 
by means of a Schechter curve and found that the characteris- 
tic dust mass at the break is M* ~ 4.7 ■ 10^ Mq at redshift z=5. 
Similar values can be achieved by models such as Lai during 
their first Gyr of evolution (we expect these galaxies to be close 
to a typical present-day L* galaxy, Fig.[T2li and are smaller than 
the ones predicted for the Hal model which, in turn, should rep- 
resent the bright end of the present-day luminosity function. 

Another class of objects which can trace the evolution 
of low- and intermediate sized ellipticals are Lyman Break 
Galaxies (LBGs). LBGs are starburst galaxies observed at high 
redshift (z > 2.5) and are identified by the colours of their 
far ultraviolet spectral energy distribution around the 912 A 
Lyman continuum discontinuity (Giavalisco 2002). The interpre- 
tation of their abundance patterns suggests that LBGs may rep- 
resent young spheroids, observed during their main SF episode 
(Matteucci & Pipino 2002, where we refer the reader interested 
in a more complete analysis), under the assumption that half 
of Fe was hidden in dust. Here we briefly repeat the study of 
Matteucci & Pipino (2002) by comparing the output of the chem- 
ical evolution code with the observed abundances in MS 1512- 
cB58, the brightest LBG known so far owing to its gravitation- 
ally lensed nature. A stellar mass of ~ IO^'M©, an age lower 
than ~ 300 Myr and a SFR of 4-QMQyr^^ (but we can consider a 
range of possible values of 20-80 Mojr"' given the uncertain- 
ties affecting the derivation of the SFR), have been reported for 
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Fig. 8. Predicted dust production rates from various sources for the chemical evolution model Lal+Itoh (fiducial case) for an 
elliptical galaxy. In the lower (upper) panel, we show the results for the carbon (silicate) dust. Solid lines: contribution by low and 
intermediate mass stars (LIMS). Dotted lines: contribution by type 11 SNe. Dashed lines: contribution by type la SNe. 




Fig. 9. Fractions locked in dust for C (solid line), O (dotted), S (dot-dashed), Fe (dashed) for fiducial models Hal (left panel) and 
Lai (right panel). 



MS1512-CB58 (Pettini et al. 2002). In order to have a realistic 
model for such a galaxy, we extend the class of models of Pipino 
et al. (2005) down to a Mium = 10 '"Mq spheroid, with a SF ef- 
ficiency V = 5 Gyr ' and an infall timescale t - 0.5 Gyr. The 
maximum value for the SFR is ~ 20^^03^"'. Since we focus on 
the very early phases of evolution, the assumptions regarding the 
dust destruction in the hot halos discussed in the previous para- 
graphs do not affect the results. At 130 Myr, our model predicts 
the following abundance ratios: [Fe/H]=-1.13, [O/H]=-0.24, as 
well as [Mg/H]=-0.62 and [Si/H]=-0.64, in remarkable agree- 
ment with the values reported by Pettini et al. (2001) and Tepliz 
et al. (2000). In particular, they found [Fe/H]=-1.10, [0/H]=- 
0.26, as well as [Mg/H],[Si/H]~ -0.42, with quoted uncertainties 
of about 0.2 dex. Pettini et al. (2002) claimed that part (roughly 



half) of the Fe could be locked into dust. On the other hand, 
the theoretical expectation (based on the S.N.) that also Mg and 
Si should be depleted, although to a lesser extent, relative to S, 
was not satisfied by their data. An answer can be given by our 
model, which, in fact, predicts that at -100 Myr the dust frac- 
tions are 50%, 37%, 32% and 7% for Fe, Mg, S and O respec- 
tively. Moreover, by means of our model we can a-posteriori 
prove the hypothesis made by Matteucci & Pipino (2002) for 
the Fe depletion was right. The only difference between their 
model and this work is in the age estimate. In particular, we sug- 
gest here an age of ~ 100 Myr, whereas the previous estimate 
was close to 30 Myr. A slightly older age can be understood in 
terms of the formation process. The model used by Matteucci & 
Pipino (2002) tightly followed the monolithic scenario, therefore 
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Fig. 10. Upper panel: Temporal behaviour of the Fe grain destruction rates for an elliptical galaxy (Pipino et al.(2005)'s model Lai). 
Dotted: Lal+MK case. Solid: Lal+Itoh case. Lower panel: Evolution of the Fe dust fraction in the ISM of the above models. 




time (Gyr) 

Fig. 11. Evolution of the Fe abundance in the hot phase of an elliptical galaxy for the Lai model. Solid lines: Lal+Itoh model; 
Dotted Line: Lai +MK model. 



the SFR was maximum at the beginning. The present model, ow- 
ing to a non-negligible infall (e-folding time of 0.5 Gyr), takes a 
longer period to assemble the same amount of stars. 



For a more general comparison with LBGs, we find that the 
value for the dust to gas ratio D - 0.090© predicted by this 
last model is in the range of values predicted by Inoue (2003). 
Moreover, concerning the quantity Dz, we estimate a value of 
0.13 (i.e. 0.2 times the S.N. value) which is consistent with the 
range O.Ol-lDz©, as predicted by Inoue (2003). At the same 
age, model Lai features Dz ~ 0.7 > Dz^q', this value might 
reflect the fact that models such as Lai begin to be too massive 
to properly represent LBGs. These results reinforce the analysis 
made by Matteucci & Pipino (2002), therefore we conclude that 
LBGs are very likely to be the progenitors of present-day low 



mass spheroids. 



3.3. Dust evolution in irregular galaxies 
3.3.1. Dust Production Rates 

In Figure [T3] we show the dust production rates for the models 
of an irregular galaxy with continuous SF and a starburst irregu- 
lar. All the dust production rates calculated for the IC model have 
a smooth behaviour and in general present a continuous increase, 
lasting several Gyrs, and then reach a plateau. On the other hand, 
the production rates calculated for the starburst irregular reflect 
the shape of the SF and present a gasping or intermittent be- 
haviour. An interesting difference between the dust production 
rates predicted for iiTegular galaxies and the one predicted for 
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Fig. 12. Upper panels: Dust to gas mass ratio (D) for Pipino et al.(2005)'s model Lai. Dotted: Lal+MK case. Solid: Lal+Itoh case. 
Lower panel: Evolution of the dust mass in the ISM for the above models. 



the S.N. model concerns the production of silicate dust by low 
and intermediate mass stars. In the case of the S.N. model, the 
Si dust production by LIMS starts immediately after the first SF 
episode. On the other hand, for the IC model Si dust production 
by LIMS starts at ~ 5 Gyr, whereas for the starburst model it 
does not start at all. This reflects the fact that in dwarf galax- 
ies and in low metallicity systems in general, it is rare to have 
Xc/Xq > 1 (see section 2.1). 



3.3.2. The dust fractions 



In Figure [141 we show the predicted evolution of the dust frac- 
tions for the IC and the starburst irregular model. In the case 
of the IC model (Fig. [141 lower panel), the dust fractions are 
maximum at the beginning and decrease progressively up to the 
present time, when C, O and S have dust fractions lower than 
~ 0.1, whereas Fe of ~ 0.3. These values are considerably lower 
than the ones predicted for the S.N. and for the present-day el- 
lipticals. In this case, the main reason is the absence of the dust 
accretion term and the continuous SF history, which gives rise to 
continuous type II SN explosions. As seen in section 3. 1 .4, these 
objects are the main responsibles for dust destruction, hence play 
a dominant role in keeping the dust fractions low in irregular 
galaxies. 

In the upper panel of Fig. [141 we show the evolution of the dust 
fractions for the starburst irregular model. This model is charac- 
terised by long periods of no SF activity. During these periods, 
all the dust fractions tend to remain constant or to increase, since 
the only contributors to destruction are type la SNe. Each star- 
burst coiTesponds to massive type II SN activity, hence to con- 
siderable dust destruction and consequent sharp decrease of all 
the dust fractions, which is visible at 1, 10 and 13 Gyr. At the 
present time, the predicted dust fraction for the starburst model 
are in general higher than the ones calculated for the IC model. 



3.3.3. Dust masses locked up in galaxies and ejected into 
the IGM 



In Figure \T5\ the thick and thin lines are the predicted evolu- 
tion of the dust masses for IC and starburst models. The solid 
lines are the dust masses locked up into the galaxies, whereas the 
dotted lines are the dust masses ejected into the IGM. Our pre- 
dictions indicate that, for an irregular galaxy of ~ 5 ■ IO^Mq, the 
present-day galactic dust masses are of the order of 1 -2- IO'^Mq, 
iiTespective of the SF history. These values are compatible with 
the dust masses observed at the present time in dwarf irregulars 
(dirr) and blue compact galaxies (BCG), which can span from 
~ lO^Mo up to ~ lO^'Mo (Lisenfeld & Ferrara 1998). Our values 
for the dust masses in local dwarfs and irregulars are compatible 
also with observational determinations achieved by means of the 
Spitzer Space Telescope. Cannon et al. (2006) derive for the lo- 
cal dwarf NGC 6822 a total dust mass of 8.3 xIO'^Mq. For dwarf 
galaxies of the M8 1 group, Walter et al. (2007) find dust masses 
of 10'* - IO^Mq, in perfect agreement with the values we predict 
for dwarfs. 



Larger diff^erences are predicted for the dust masses ejected 
by IC and starburst irregular galaxies into the IGM. Our results 
indicate that, for the IC model, the galactic wind develops very 
early, at ages comparable to ~ 1 Gyr, whereas for the starbust 
irregular the onset of the wind occurs after 10 Gyr of evolution. 
IC and starburst irregulars eject into the IGM ~ 2 • IQ^Mq and 
~ 1 • IO'^Mq of dust, respectively. Calura & Matteucci (2006a) 
have shown that the bulk of the IGM metal enrichment is due to 
giant spheroids, despite the fact that the number density of dwarf 
irregulars is larger than the one of the giant ellipticals. Since 
the ejected dust mass is proportional to the total ejected metal 
mass, we can conclude that ellipticals should be considered as 
the major contributors to the dust enrichment of the IGM. 
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Fig. 13. Predicted dust production rates from various sources for a chemical evolution model of an irregular galaxy with continuous 
SF (left) and a starburst irregular galaxy (right). In the lower (upper) panels, we show the results for the carbon (silicate) dust. Solid 
lines: contribution by low and intermediate mass stars (LIMS). Dotted lines: contribution by type II SNe. Dashed lines: contribution 
by type la SNe. 




Fig. 14. Time evolution of the fractions in dust for various elements for the IC model (lower panel) and for the starburst model 
(upper panel). The Solid lines, dotted lines, dashed lines and dot-dashed lines represent the evolution of the C, O, Fe and S dust 
fractions, respectively. The evolution of the dust fractions for Mg and Si is identical to the Fe one. 



3.3.4. Dust to gas and Dust to metals ratios 

In the case of irregular systems, in general the observed dust to 
gas ratio is defined as: 



Mhi 

where M^/„j, and Mui are the dust and HI mass determined 
from the far infrared and from the radio emission, respectively 
(Lisenfeld & FeiTara 1998). In the models, the neutral H mass is 
Mhi - Mgas ■ Xh, where M^yas is the total gas mass and Xh is the 
H mass fraction. In Figure □11 we show the D' calculated for the 
IC model and for the starburst model as a function of time (upper 
left panel) and as a function of the metallicity (lower left panel). 



and compared to a set of observations in local dwarf galaxies. At 
any time, the IC model is characterised by higher D' values than 
the starburst model. 

In the lower panel of Figure [161 a proxy of the metallicity 
we use the quantity 12+log(0/H). The observational values used 
here are taken from a compilation by Lisenfeld & Ferrara (1998) 
(see caption of Fig. [T6] l, and refer to BCD (solid circles) and 
dliT (open squares) galaxies. The dust to gas ratio of dwarf ir- 
regulars and BCDs as a function of the metallicity was studied 
by Hirashita (1999) and Hirashita, Tajiri & Kamaya (2002). In 
particular, Hirashita et al. (2002) suggested that an intermittent 
SF history allows one to reproduce the observed dispersion in 
the relation between the dust to gas ratio and the metallicity. 
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Fig. 15. Time evolution of the dust masses for the IC (thick lines) and starburst irregular (thin lines) models. Solid Unes: dust mass 
locked up into the galaxies. Dotted lines: the dust masses ejected into the IGM. 
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Fig. 16. Left: Evolution of the dust to gas ratio calculated for the IC model (solid lines) and for the starburst model (dashed lines) 
as a function of time (upper panel) and metallicity (lower panel). In the lower panel, the solid circles and open squares are D' values 
observed in local dliT and BCD galaxies, respectively (Lisenfeld & FeiTara 1998). The solid and empty inverted triangles represent 
upper limits to the observations for local dlrr and BCD galaxies, respectively, calculated by assuming an uncertainty of a factor of 
4. Right: Evolution of the dust to metals ratio calculated for the IC model (solid lines) and for the starburst model (dashed lines) as a 
function of time (upper panel) and metallicity (lower panel). The thick and thin lines represent the predicted D' calculated assuming 
for the condensation efficiencies the values suggested by D98 and a constant value of 0.1, respectively. All the symbols are as in the 
left lower panel. 



The starburst model has a present-day metallicity lower than the 
IC model. We note that the predicted present-day values are in 
general higher than the majority of the observed values. The 
compilation of observations presented by Lisenfeld & FeiTara 
(1998) is based on IRAS data. The observed determinations of 
the dust mass are likely to represent lower limits to the actual val- 
ues, owing to the undetected cold dust component (D98, Popescu 
et al. 2002). Lisenfeld & Ferrara (1998) estimate a total error in 
the dust to gas ratio to be of a factor of 4, taking into account all 
the possible factors of uncertainty (contribution from very small 



grains, cold dust, molecular gas and variations in the HI and 
optical diameter). By taking into account this uncertainty (solid 
and open inverted triangles in Fig.fTSTi, the observed dust to gas 
ratios are consistent with our predictions for both irregular types. 
The dust to metals ratio D'^ observed in dwarf irregulars can be 
defined as: 



Mdu 



HI • ■ 



D' 

Y 



(20) 



In section 3.1.2, we have seen that the dust depletion patterns 
(and hence also the dust to metal ratio) are determined by the 
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balance between destruction and accretion and that, in the S.N., 
it is impossible to derive constraints on the dust condensation 
efficiencies from the analysis of the observed depletion pattern. 
In dwarf irregulars, dust accretion is likely to play a negligible 
role or to be absent, hence the variation of the condensation 
efficiencies has important effects on the dust to metal ratio. In 
the right upper and lower panels of Figure [16] we show the pre- 
dicted evolution of D'^ in dwarf irregular and starburst galaxies 
as a function of time and metallicity, respectively, calculated 
assuming the set of condensation efficiencies as suggested by 
D98 (thick Hnes) and a constant value of 6^.^ = d'" = S\' = 0.1 
(thin lines). At variance with the S.N., the adoption of different 
condensation efficiencies has noticeable effects on the dust to 
metal ratio. In particular, the D'^ values calculated by assuming 
the D98 condensation efficiencies are larger than the ones 
predicted by assuming a constant value of 0. 1 . 

In the lower right panel of Fig. [16] we show the predicted 
evolution of D'^ as a function of the metallicity for the IC model 
(solid lines) and for the starburst model (dashed lines). The 
big open circles and the big open squares are the predicted 
present-day values for D'^ for the starburst and IC model, respec- 
tively. By comparing the different predictions for D^, calculated 
assuming two different sets of dust condensation efficiencies, we 
can note that in principle, the measure of the dust to metal ratio 
and of the dust depletion pattern in dwarf irregulars could allow 
one to put solid constraints on the dust condensation efficiencies. 



4. Conclusions 

By means of chemical evolution models for galaxies of diff'er- 
ent morphological types, we have performed a study of the cos- 
mic evolution of the dust properties in different environments. 
By adopting the same formalism as developed by D98, we have 
studied the evolution of the dust content of the S.N., confirm- 
ing the main results of D98. We have carried on a deep study 
of the space of the parameters used to model dust evolution and, 
thanks to the uptodate observations available in the solar vicinity, 
performed a fine tuning of the parameters. In the following, we 
outline only the new aspects which have been investigated here 
and which are different than the ones considered by D98. We 
have extended our study to ellipticals and dwarf irregular galax- 
ies, for which dust evolution has been calculated by means of 
chemical evolution models relaxing the instantaneous recycling 
approximation, namely taking into account the stellar lifetimes. 
The main results concerning the present-day dust properties for 
the galaxies studied here are presented in table 3. Our main re- 
sults can be summarised as follows. 

1) For each chemical element, we have compared the pre- 
dicted present-day fractions to the ones observed in the Local 
Interstellar Cloud. This analysis has been useful to test the vari- 
ous parameters involved in our study, i.e. the dust condensation 
efficiencies 5?'*', 6\", 6'' and the dust destruction efficiency e. In 
the most realistic physical situation, i.e. once both the dust de- 
struction and accretion processes are taken into account, for all 
the elements the dust fractions are nearly independent from the 
choice of the dust condensation efficiencies. The main processes 
determining the gas fractions are dust accretion and dust destruc- 
tion. We can reproduce the observed dust fractions by assuming 
that the dust destruction efficiency depends on the properties of a 
given chemical element. A physical justification for this assump- 
tion might be found in the different condensation temperatures 
Tc of the elements. These assumptions allow us to reproduce also 



the gas and cosmic abundances observed in the Local Interstellar 
Cloud. 

2) For all the elements, the evolution of the dust fractions is 
strongly determined by the SF history. At the beginning of each 
SF event, the dust fractions present a sharp peak, immediately 
followed by a very steep decrease, due to dust destruction by in- 
tense type II SN explosions. 

During most of the cosmic time, i.e. from ~ 3 up to 10 Gyr, in 
the S.N. the dust fractions rise monotonically. The effect of the 
SF threshold causes rapid oscillations in the dust fractions only 
in the last 3 Gyr of evolution. 

3) For the dust destruction rate, we predict present-day values 
of ~ 5 X 10"''yr"' and ~ 8 x 10"''yr"' for silicate and carbon 
dust, respectively. Similar values are predicted for the accretion 
rates for both dust types. Our values are in good agreement with 
the estimates by Draine (1990) and Tielens (1998) for a cold 
medium with different approaches. 

4) For both type la and type II SNe, we compare the production 
and destruction rates and how these quantities vary throughout 
the evolution of the S.N. For type II SNe, dust destruction domi- 
nates over dust production throughout almost all the cosmic his- 
tory. On the other hand, for type la SNe during the first Gyr of 
galactic evolution. Si dust production has dominated over dust 
destruction, but with negligible consequences on the total Si dust 
mass. 

5) We study the evolution of the dust to gas and dust to met- 
als ratio as a function of time and metallicity in the S.N. and 
in the outermost regions of the Galactic disc. The disc outskirts 
evolve with a SF activity less intense than the S.N., with a lower 
present-day dust to gas ratio and a higher present-day dust to 
metals ratio. For both quantities, we compare the value predicted 
by means of the S.N. model to the one observed in the Galaxy, 
finding a very good agreement. 

6) In elliptical galaxies, type la SNe are the major dust factories 
in the last 10 Gyr. With our models, we successfully reproduce 
the dust masses observed in local ellipticals (~ lO^M©) by means 
of recent FIR and SCUBA observations. 

7) We have shown that, by considering a reduced dust destruc- 
tion in a hot and rarefied medium, dust is helpful in alleviating 
the iron discrepancy observed in the hot gaseous halos surround- 
ing local ellipticals (Arimoto et al. 1997). In this medium, we 
predict a Fe abundance of [Fe/H] =0.1 5-0.6 (depending on the 
adopted prescription for dust destruction) whereas without dust 
Pipino et al. (2005) found [Fe/H]> 0.85. Moreover, the inclu- 
sion of dust improves the agreement between the predicted X- 
ray temperature (0.7 keV) and luminosity (~ 10''^ergs"') for a 
model of a giant elliptical and the observations. 

8) The dust masses observed in high-redshift SCUBA galax- 
ies, the most likely progenitors of the local giant ellipticals, 
are successfully reproduced by our models, which predict up to 
~ IQ^^'^Mq of dust during the high-redshift star forming phase. 
Furthermore, the dust treatment is very helpful in reproducing 
the abundances and the dust to gas ratios observed in Lyman 
Break Galaxies. In the specific case of the galaxy MS1512-cB58, 
we confirm the assumptions on dust made by Matteucci & Pipino 
(2002) and we reinforce their conclusions, i.e. that LBGs are 
very likely to be the progenitors of the local low mass spheroids. 

9) As for the metals, we predict that ellipticals play a major role 
in the dust enrichment of the IGM/ICM during galactic winds. 
A minor fraction (i.e. ~ 10""* of the dust mass ejected by ellipti- 
cals) is ejected by Irregulars. 

10) The two models used to study dust evolution in dwarf irreg- 
ular galaxies present very different dust production rates. For the 
IC model (irregular galaxy with continuous star formation), the 
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Table 3. 

Predicted present-day dust properties for the galactic models used in this work. The various models are listed in column 1. In columns 2, 3, 4, 5, 6 
we present the predicted present day dust production rates, dust destruction rates, dust accretion rates, dust masses, dust to gas ratios and dust to 

metals ratios, respectively. 



production rates have a smooth behaviour throughout the whole 
evolution of the system. For the starburst model, the evolution 
of the production rates reflects its intermittent SF history. Owing 
to the different SF histories, the starburst model presents higher 
dust fractions than the IC model at any time. 

11) The predicted present-day dust to gas ratios for irregular 
galaxies are in general higher than the majority of the observed 
values. However, the determination of the dust mass in these 
galaxies is affected by several sources of uncertainty. By assum- 
ing a factor of 4 uncertainty in the observed data, as suggested 
by Lisenfeld & Ferrara (1998), the observed dust to gas ratios 
are consistent with our predictions. 

12) In dwarf irregulars, dust accretion is likely to play a negligi- 
ble role or to be absent, hence the variation of the condensation 
efficiencies has important effects on the dust to metal ratio. In 
principle, a precise determination of the dust to metal ratios and 
of the depletion pattern in dwarf irregular galaxies could be help- 
ful to put solid constraints on the dust condensation efficiencies. 

13) Finally, it is interesting to compare our results for the 
present-day dust masses predicted for elUpticals and irregulars 
with the recent observational determinations, possible thanks to 
the Spitzer Space Telescope. 

In local ellipticals, Kaneda et al. (2007) and Panuzzo et al. 
(2007) derive dust masses typically of the order of 10^ - IO^Mq. 
The values computed by means of our fiducial model are com- 
patible with the observationally derived values quoted above (see 
table 3). 

Also our predictions for the dust masses of dwarfs/irregular 
galaxies are in agreement with the Spitzer results. In fact, for 
local dwarf galaxies, dust masses of lO"* - IO^Mq have been de- 
tected (Cannon et al. 2006, Walter et al. 2007), in perfect agree- 
ment with the values we predict for present day dwarf galaxies, 
i.e. 1 - 2 X IO^Mq. 
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